Anomalous diffusion and stretched exponentials in heterogeneous glass-forming 

liquids: Low-temperature behavior 
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We propose a model of a heterogeneous glass forming liquid and compute the low-temperature 
behavior of a tagged molecule moving within it. This model exhibits stretched-exponential decay 
of the wavenumber-dependent, self intermediate scattering function in the limit of long times. At 
temperatures close to the glass transition, where the heterogeneities are much larger in extent than 
the molecular spacing, the time dependence of the scattering function crosses over from stretched- 
exponential decay with an index b = 1/2 at large wave numbers to normal, diffusive behavior with 
& = 1 at small wavenumbers. There is a clear separation between early-stage, cage-breaking 13 
relaxation and late-stage a relaxation. The spatial representation of the scattering function exhibits 
an anomalously broad exponential (non-Gaussian) tail for sufficiently large values of the molecular 
displacement at all finite times. 



I. INTRODUCTION 

A growing body of evidence, both experimental 
and numerical 0, |^ , points to the existence of intrinsic 
spatial heterogeneities in glass forming liquids at tem- 
peratures slightly above the glass transition, and the rel- 
evance of these heterogeneities to stretched exponential 
decay of correlations and anomalous, non-Gaussian diffu- 
sion. Here we propose a simple model of a heterogeneous 
glass former and use it to compute observable proper- 
ties of a tagged molecule moving within it. In particular, 
we show that this model naturally predicts stretched- 
exponential decay when the scale of the heterogeneity is 
much larger than the molecular spacing; and we compute 
deviations from Gaussian displacement distributions dur- 
ing intermediate stages of the relaxation process. 

Our model is loosely motivated by the excitation-chain 
(XC) theory of the glass transition proposed recently by 
one of us (JSL) 0, H This theory suggests that a glass- 
forming liquid at a temperature T not too far above the 
glass transition temperature Tq consists of fluctuating do- 
mains of linear size R*(T), within which the molecules 
are frozen in a glassy state where they have little or no 
mobility. R*{T) is a length scale that characterizes dy- 
namic fluctuations in the equilibrium states of these sys- 
tems. The theory predicts that R*{T) diverges near Tq 
like (T — To)^^ and decreases to zero at the upper limit 
of the super- Arrhenius region, T — Ta- 

Characterizing the more mobile material that lies out- 
side the hypothetical frozen domains remains one of the 
deeper challenges in glass physics. In the language intro- 
duced by Widmer-Cooper and Harrowell [9|, the mobile 
molecules reside in regions of high "propensity." Those 
authors perform molecular dynamics simulations of a 
glasslike, two-dimensional, binary mixture. They average 
the mean-square displacement of each molecule over an 
isoconfigurational ensemble in which the initial positions 
of all the molecules are fixed, but the initial velocities are 
selected at random from a Boltzmann distribution. The 
propensity of a molecule is its mean-square displacement 



during a structural relaxation time. Contour maps of 
propensity as a function of initial molecular positions do 
indeed exhibit well defined domains of low and high mo- 
bility. 

Widmer-Cooper and Harrowell compare their propen- 
sity maps with corresponding plots of the potential en- 
ergy, which is taken to be indicative of the local struc- 
ture. Interestingly, the propensity appears to be almost 
entirely uncorrelated with the potential energy. In retro- 
spect, this lack of correlation between short-range struc- 
ture and mobility may not be surprising. Such a corre- 
lation would be expected for a slowly coarsening system 
consisting of domains of two structurally distinct phases, 
for example, a slightly supersaturated (mobile) fluid in 
which (frozen) crystallites are growing. In that case, the 
domains are very nearly but not quite in thermodynamic 
coexistence with each other. A true glass-forming liquid, 
however, is in a state of thermodynamic equilibrium. The 
temperatures and chemical potentials are uniform every- 
where, independent of propensity; and the local pair cor- 
relations, which ought to be determined primarily by the 
temperature in an equilibrium state, appear to be statis- 
tically uniform as well. 

A local quantity that does correlate with the propen- 
sity is the mean-square vibration amplitude of a molecule 
averaged over times much longer than its period of oscil- 
lation but much shorter than the time required for irre- 
versible rearrangements to occur. A large mean-square 
vibration amplitude - or, equivalently, a small Debye- 
Waller factor - implies that a molecule is participat- 
ing in a soft, low frequency, elastic mode. It is well 
known that soft modes are abundant in amorphous mate- 
rials near jamming transitions. jlflj Moreover, it appears 
(for example, in Fig. 4 of [l^) that the largest displace- 
ments in such modes lie predominantly along one dimen- 
sional, chainlike paths. Thus, regions of high propensity 
may be elastically soft, and may somehow be correlated 
with excitation-chain activity. Of course, local variations 
of elastic stiffness are ultimately structural in nature; 
but soft modes are collective phenomena involving many 
molecules, and are not easily detected by measuring near- 
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neighbor pair correlations. 

In the model proposed here, glassy domains of low 
propensity are separated from each other by regions of 
higher propensity in which many of the molecules are 
mobile. This domain structure undergoes persistent fluc- 
tuations on time scales of the order of r* = R*'^ /Da^ 
where Da = P'jTa is the diffusion constant associ- 
ated with a relaxation, I is approximately the average 
molecular spacing, and Tq, is the strongly temperature- 
dependent, super- Arrhenius, a relaxation time. In con- 
trast to the slow fluctuations of domain boundaries, indi- 
vidual molecules within high-propensity regions are ran- 
dom walkers with a diffusion constant Dm = (-"^ /tm, 
where tm is a temperature-dependent time scale ~ pre- 
sumably not super- Arrhenius - that characterizes mobile 
molecular displacements. 

To a first approximation, a tagged molecule at any 
given time is either frozen in a glassy domain or is mobile. 
If the former, then it waits a time of the order of r* before 
it is encountered by mobile neighbors and itself starts 
to undergo displacements. If the latter, then it diffuses 
for some distance until it finds itself frozen again in a 
glassy domain. That distance is at least a few molecular 
spacings t, but it seems more likely to be proportional to 
R* if the domains are large and if diffusive hopping takes 
place more readily within the mobile regions than across 
their boundaries and into the glassy domains. Therefore, 
we assume that the time elapsed during mobile motion 
is of the order of R*^/Dm = {R*lifTM- The ratio of 
these two time scales is t^Itm = A. At temperatures 
well below Ta, we expect that A ^ 1. 

A similar picture of diffusion in fluctuating regions of 
varying mobility emerges in kinetically constrained mod- 
els [n], and in a recent model of a gel [l3|. Our 
analysis is most closely related to the diffusion model of 
Chaudhuri et al which has been a valuable start- 
ing point for the present investigation. We go beyond 
14| by including a mechanism for producing late-stage 
stretched-exponential relaxation, and by distinguishing 
that mechanism from the diffusion of mobile molecules. 
In contrast to ||14], our molecules switch back and forth 
between frozen and mobile behavior and thus explore the 
geometry of the domain structure. As a result, our model 
exhibits a clear separation between the slow a relaxation 
and the faster mobile motions that we interpret as /? re- 
laxation. It also respects time-translational symmetry, 
which seems to be violated in [l^ . 

The preceding discussion pertains just to the case of 
large-scale heterogeneities at temperatures only a little 
above To and well below Ta- Some of the most interesting 
and experimentally accessible physics, however, occurs 
near the crossover to small-scale heterogeneity near Ta- 
Constructing a theory of diffusion in the crossover region 
requires a temperature-dependent analysis, which will be 
described in a following paper. 



II. BASIC INGREDIENTS OF THE MODEL IN 
THE LIMIT OF LARGE-SCALE 
HETEROGENEITY 



As outlined above, there are two different mechanisms 
to be described probabilistically in this sytem, each op- 
erating on its own characteristic time scale. The separa- 
tion between time scales, characterized by a large value 
of A, suggests that the natural mathematical language in 
which to discuss this system is that of a continuous time 
random walk (CTRW) in which a tagged molecule al- 
ternates between long waits in glassy domains and faster 
displacements in mobile regions. In this Section, we com- 
pute the waiting time distributions for molecules in glassy 
and mobile domains, and the corresponding probability 
distributions for the mobile displacements. 

We immediately encounter a problem, however, be- 
cause there are more than just two characteristic time 
scales in this system. Analog experiments [3] and numer- 
ical simulations [1] typically follow the motions of parti- 
cles starting from their initial positions, and resolve these 
trajectories on length scales smaller than the interparti- 
cle spacing, and on time scales of the order of the shortest 
vibrational periods. This is a continuous range of scales. 
At its long-time limit, it includes the slow modes that 
are involved in molecular rearrangements, whose ampli- 
tudes may exceed the Lindemann melting criterion in the 
mobile regions. Thus, the cage-breaking, /^-relaxation 
mechanism may be part of the same continuous range of 
time scales that includes the intra-cage vibrational mo- 
tions. Indeed, no sharp distinction between vibrational 
and cage breaking time scales seems to appear in exper- 
imental data. If these are not clearly distinct kinds of 
events, then the CTRW approximation is not valid in 
the short-time limit. 

In view of this difficulty, we choose here to consider 
only the (5 and a time scales, tm and Ta respectively. By 
neglecting the shorter time scales, i.e. those of the order 
of tq ^ femtoseconds - picoseconds, we imply that we are 
not resolving length scales much smaller than the molec- 
ular spacing i. We then assume that the probability that 
a molecule has left its cage at times of order tm ^ tq is 
determined by the probability that it has mobile neigh- 
bors and participates in local molecular rearrangements. 
At the longest time scales Ta ~ seconds (by definition, 
the viscous relaxation time at the glass temperature), 
the probability that a molecule has escaped from a large 
glassy domain is determined by the probability that it 
has been encountered by the boundary of the domain. 

Consider first the slowest motions, i.e. those associated 
with fluctuations of the domain boundaries on length 
scales R* and time scales t*. Define t* = t/T^, where 
t is the physical time in seconds; and let ipait*) be the 
normalized probability distribution for the time that a 
molecule spends in a glassy domain before entering a mo- 
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bile region. Write this distribution in the form 

^ain^ / dpWip) (2.1) 

Jo P 

where p is the hnear size of a domain in units R* . W{p) 
is a normalized distribution over these sizes, and the re- 
maining factor inside the integrand is a normalized dis- 
tribution over t* . The quantity appearing in the 
exponential in Ea. (|2.ip is the lowest eigenvalue of the 
diffusion kernel for a molecule moving in a domain of 
size p. Note the similarity to the trapping models dis- 
cussed, for example, in [l^, [11, |T3|; but also note that 
the "trap" here is a two dimensional subspace bounding 
a three dimensional domain. The quantity that is fluc- 
tuating in a normal diffusive manner is the mean square 
of the distance between the tagged molecule and the do- 
main boundary. Because the boundary surrounds the 
molecule, the direction of the diffusive drift is irrelevant; 
any point of contact on the boundary is equivalent to any 
other. Therefore, it makes little difference whether the 
molecule is the diffuser and the boundary is the target, or 
- as in this case - the boundary diffuses and the molecule 
is the target. 

Strictly speaking, Eq. (j2.ip is a long-time approxima- 
tion; at shorter times, the higher eigenmodes of the diffu- 
sion kernel make non-negligible contributions. However, 
this approximation is qualitatively adequate for our pur- 
poses at shorter times as well because ipcit*) is well be- 
haved at small t* and, as mentioned previously, we are 
not trying to include the short-time behavior in this func- 
tion. If we assume that the distribution over values of p 
is Gaussian, W{p) oc p^ exp(— p^) in three dimensions, 
then 

^PG{t*)^2e-^^. (2.2) 

That is, we find a rudimentary but nontrivial stretched 
exponential of the form exp (—const, x t*'') with b — 1/2. 

There is little reason to expect that W{p) remains 
Gaussian out to the large values of p that determine 
the long-time behavior of 4'G{t*)- Deviations from the 
Gaussian would produce different indices. For example, 
as the temperature increases and the glassy domains be- 
come smaller, their size distributions might cut off more 
sharply than a Gaussian. Thus, if W{p) oc exp(— p™), 
with m > 2, then h = m/ {m + 2) — > 1 as m ^ oo. The 
problem of computing this distribution or, equivalently, 
values of m from first principles may eventually become 
solvable as we learn more about the statistical mechan- 
ics of glass forming liquids; but that problem is beyond 
the scope of this investigation. For present purposes, the 
important points are that a plausible distribution W{p) 
produces stretched-exponential behavior, and that the 
mechanism by which this happens could produce a tem- 
perature dependent index h. 

Next, consider the waiting-time distribution ^pM{t*) 
for a molecule in a mobile region. This situation is qual- 
itatively different from that of a molecule in a glassy do- 
main because the distance traveled by the molecule is 



related to the time during which it remains mobile. As 
stated earlier, we assume that a mobile molecule diffuses 
a distance of order R* before reentering a glassy region, 
so that its residence time in the mobile region is of order 
tm {R* /£)'^ which, in t* units, is simply A~^. Then, for 
simplicity, assume an exponential waiting-time distribu- 
tion: 

V'Af(t*) = Ae-^**. (2.3) 

Compared to ipG{t*) in Eq. (|2.2p . ipM{t*) is sharply 
peaked near t* = if, as expected, A is large. To com- 
plete the model of mobile motion, we need the conditional 
probability pMir* , t*) for diffusion over a scaled distance 
r* = r/R* in time t* . The natural choice is 

PMir\n . (^^^^ exp (-^) , (2.4) 

which is a normalized, three dimensional distribution 
over r* . 

III. CONTINUOUS-TIME RANDOM WALKS 

The next step is to translate these physical ingredients 
of the model into the language of continuous-time ran- 
dom walks. p^[20j Define two different probability distri- 
bution functions, nalr* , t*) and 71m(?'*, t*), for molecules 
starting, respectively, in glassy domains or mobile re- 
gions, and moving distances r* in times t*. Each molec- 
ular trajectory consists of a sequence of transitions be- 
tween glassy domains and mobile regions. A single tran- 
sition in which a molecule starts in a glassy domain at 
time i* and ends in a mobile region at time ^2 (without 
having changed its actual position) occurs with probabil- 
ity TpG (^2 ) • Similarly, a transition in which a molecule 
starts in a mobile region at time ti and position r*, and 
ends in a glassy domain at time ^2 s-nd position rj, has 
probability -^a/ (tj - i * ) pi/ ( | r2 - r J | , t2 - t ^ ) ■ The prob- 
ability that a molecule starts or arrives in a glassy do- 
main at time tl and is still there at the final time t* is 
4>Git* — tl), where 

/•oo 

Mn^J^ M*')dt*' = {1 + 2 VF)e~^^. (3.1) 

Finally, the probability that a molecule starts or arrives 
in a mobile region at time and position and is still 
in that region at the final time t* and at position r*, is 
<l>M(t* - tl)pM{\r* - rl\,t* - tl), where, in analogy to 
Eq.dSIl), 

<l>M{t*)^J^ iJM{tl)dtl=e-^''. (3.2) 

The multiple convolutions that describe each trajec- 
tory are, as usual, converted into products by computing 
Fourier-Laplace transforms of each function: 

nj{k,u) = J e-'^-^'dr* e-""'' rijir* ,t*) dt* , (3.3) 
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where j = G, M. For jumps starting in glassy regions, 
we need: 



^g{u) = I e 

'0 



-ut* 



and 



(f)c{u) 



1 - -ijjGiu) 



(3.4) 



(3.5) 



For the mobile regions, the spatial Fourier transform of 
the conditional probability distribution pm{t* ^t*) is 



PMik,t*) = exp - 



(3.6) 



therefore, define 

fM{k,u) E 



e-"**PM(fc,r)VAf(r)dr 



A 



and 



gM{k,u) 



M + A(l + fc2/2)' 

e-"**PM(fc,t*)(/'M(t*)rft 
1 



(3.7) 



u + A(l + A:72) 



(3.8) 



Putting these pieces together, and summing over in- 
definitely many individual jumps in each trajectory, we 
find 



nG(k,u) 



and 



nM{k,u) 



(j)G{u) + ffAf (fc, u) rpoju) 

1 - ■0g(m) fM{k,u) 
1 NG{k,u) ^ 
u W{k, u) 

4>g{u) fM{k,u) + gM(k,u) 

1 - ■0g(w) fM{k,u) 
1 NM{k,u) 



(3.9) 



u W{k,u) ' 

where 

Naik, u)=[l- i^aiu)] (1 + fcV2) + u/A; 



and 



A^m(A:,u) = 1-Vg(m)+m/A; 



W{k, u) = l- ^g{u) + k^/2 + m/A. 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



In both Eqs. (|3.9p and (|3.10p . in the numerators of the 
first expressions, the first term corresponds to trajecto- 
ries that end in glassy domains, and the second to those 
that end in mobile regions. 




FIG. 1: Intermediate scattering function Fs{k,t*) with Vc = 
0.5 and A = 100, for k = 0.8, 1, 1.2, 1.5, 2, 3, 5, 10, reading 
from top to bottom. 



IV. INTERMEDIATE SCATTERING 
FUNCTION 

The conventional way to study the diffusion mecha- 
nisms discussed here is to measure the self intermediate 
scattering function Fs{k,t), which, in the present no- 
tation, is the mixed Fourier-time representation of the 
weighted average of no and um ■ That is, 

Fsik,t*)^rGnGik,t*) + {l-rG)nMik,t*), (4.1) 

where Vg is the probability that a molecule starts its 
motion in a frozen, glassy environment; and 



nj{k, t*) 
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^e"** nj(fc,M), j = G,M. (4.2) 



The inversions of Laplace-transforms in Eq. (l4.2p are 
nontrivial because the stretched-exponcntial function 
^-)g{u), defined in Eq. (j3.4p . has an essential singularity 
at u = 0. For t* > 0, the integrations over u in Ea. (|4.2p 
must be performed by closing the contour in the negative 
half It-plane, which requires that hj{k, u) be analytically 
continued to points u = |m| exp(i0). To evaluate ipdu) 
at such points, using Eq. (|3.4p . let t* = and then ro- 
tate the y contour away from the positive real axis to 
a line y = (, exp{~i9/2), < ^ < oo so that the inte- 
grand always decreases rapidly at infinity. For example, 
for u —w zLie, 9 = ±7r. 



il^Gi—w ± ie) 



-4 / id^e 



1,3/2 



T ^ 
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A{w) ^iB{w) 



(4.3) 



where Erfi is the imaginary error function. The imagi- 
nary part of ipGi~w ± ie) is the discontinuity across a 
cut along the negative u axis. With this formula, it is 
straightforward to compute the discontinuity across the 
cut for the hj{k,u) in Eqs. (|3.9[) and (|3.10[) . and in this 
way to compute the scattering functions by closing the 
contour around this cut. That is: 



hj{k, t*) 



dw 

'^w 



X Im 



Nj{k,-w + iO) 



W{k,-w + iO) 



J = G, M. (4.4) 



The integrands are well behaved as w — > 0+, therefore it 
is possible to let the lower limit of integration he w = 0. 

The graphs in Fig[l] show Fs{k,t*), computed nu- 
merically from Ea. (|4.4p . as a function of log]^g(<*), for 
A = 100, Vg = 0.5, and for a sequence of wavenumbers 
k. These are low-temperature scattering functions, where 
A 3> 1, which means that molecules spend much longer 
times frozen in glassy domains than they do when moving 
in mobile regions. The choice Vg = 0.5 is made primar- 
ily for clarity; but it seems likely that this value of Vg 
describes a system that is well into its non-Arrhenius, 
anomalous-diffusion regime. The set of curves in FigfT] 
closely corresponds to those shown, for example, in 
Howewver, as explained in Sec|TTl we do not claim to 
resolve the small-time behavior accurately. 

For the larger values of k shown in FigdJ the scat- 
tering functions exhibit two-stage relaxation. The first 
drop, which we identify as the P relaxation, occurs be- 
cause the initially mobile molecules are diffusing beyond 
their cages. Then, after much longer times indicated by 
the plateaus, the scattering curves cross over to typical a 
behavior in which all molecules, independent of whether 
or not they were initially mobile, make slow transitions 
back and forth between mobile and frozen states. This 
two-step relaxation disappears at small k, where the scat- 
tering functions are averages over distances larger than 
R* and exhibit normal, diffusive behavior. 

The most interesting feature of these scattering func- 
tions is that they exhibit a continuous range of stretched 
exponential relaxation modes. To see this, we next look 
at the long-time behavior of fiGik, t*) for relatively large 
values oit* in the range 1 < t* < 1000, and for a sequence 
of values of k. In FigO we plot — log^Q [— log^o naik, t*)] 
as a function of \ogiQ{t*), so that the slopes of the curves 
are equal to (minus) the stretched-exponential index b. 
The results are again for A = 100. 

Each of the curves in Figl2] has a constant slope over 
about two decades in indicating apparently well de- 
fined values of binfiGoa exp {~t*^). In the hmit of large 
k, nG{k,t*) is indistinguishable from the glassy waiting 
time distribution (j)G{t*) defined in Eq. ()3.ip and shown 
by the red curve in the figure; and b is accurately equal to 
1/2 for large t* . b increases toward unity as k decreases 



toward values of order unity or smaller. For any nonzero 
fc, these curves also exhibit a crossover from 6 > 1/2 at 
smah to & = 1/2 at large t* . 




FIG. 2: (Color online) Intermediate scattering function 
hcikjt*) for k = 0.7, 1, 1.5, 5, 10, reading from top to bot- 
tom. The red curve is the glassy waiting time distribution 
ipcit"), which in this plot is indistinguishable from the large- 
k limit of nG{k,t*). 

To sec this behavior analytically, we can deduce from 
Ea. (|4.3p that A{w) « 1 -t- 3?i;/2 for the very small values 
of w that are relevant at very large t* . Therefore, for 
large A, 

Vl^(fc,-u; + i0) w y - ^ +iB(w); (4.5) 

and, because B{w) <C 1 for w <C 1, the integrand in 
Ea. (|4.4p is sharply peaked at w = fc^/3 so long as fc^ <C 1. 
For t* ^ 1, the integrand in Eq. (|4.4p has another sharp 
peak at the maximum of the function exp {—wt*) B{'w), 
i.e. a.iw = w*^ l/\/t*. If fcVS < l/Vt* < 1, then the 
diffusive peak is dominant, and fiG ^ exp (— fc^ t*/3). On 
the other hand, if \/^/t* < k^/3 <^ 1, then the anoma- 
lous peak at w* is dominant, and a saddle-point esti- 
mate yields ug ^ 2 \/t* exp {—2 y/t*) as expected from 
Eq. (|3.ip . Thus, at any fixed large time t*, ug becomes 
purely diffusive in the limit of small k. Conversely, at 
fixed small k, the diffusion becomes anomalous at large 
enough t*. 

A complementary set of behaviors is illustrated 
in FiglJl which is the analog of Figl^] for the 
initially mobile molecules. That is. Fig [3] shows 
— logiQ [— logio "-M(fc, t*)] as a function of logiQ(t*), but 
for a broader range of times extending in two cases 
down to t* ~ 10^^. At large times, the preceding 
analysis remains valid; for large fc, hm ^ (t'cit*) ~ 
2 Vt* exp(-2 Vt*). As may be expected, however, and 
consistent with the behavior seen in Fig[l] um decays 
diffusively at small times (with slope —1 in this graph). 
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log^g(t*) Displacement x* 



FIG. 3: (Color online) Intermediate scattering function 
nM(fc,t*) for k — 0.7, 1, 2, 5, reading from top to bottom. 
The red curve is the glassy waiting time distribution 4>G{t*)- 



and then enters a flat plateau before crossing over to the 
long-time anomalous behavior characteristic of both of 
these scattering functions. For large enough k and/or A, 
the early-time diffusive behavior can be deduced analyti- 
cally. In this case, the relevant values of w are large, and 
A{w) « —2/w becomes negligibly small. Therefore, 

W{k,-w + iQ)Kl + ^ ~^ + iB{w), (4.6) 

where, for large w, B{w) « 2 yjl: /w'^^'^ again becomes 
small. Now the integrand has a sharp peak at w — A (1 + 
/c^/2). This peak dominates the integrand at large w for 
fiM, but its amplitude is smaller by a factor A~^ for 
fiG- Thus, for small t* , um ^ exp[-(l -I- fc2/2) At*] is 
normally diffusive, while no such behavior occurs in fiQ. 



FIG. 4: Weighted spatial distribution function Fs(x* ,t*) as a 
function of the scaled displacement x* for scaled times t* — 
0.03, 0.1, 0.3, 1.0, 3.0, 10.0, reading from left to right along 
the bottom of the graph. 



Similarly, 



nM{r*,n = -J—jLTM{r*X), (5.4) 
4 7r7'* or* 



where 



TMir*,t*) 



+ ZOO 



du 
27r iu 



k{u) e 



— k{u) r* 



(5.5) 



It is convenient to project these distributions onto, say, 
the X* axis, that is, to integrate out the perpendicular 
directions. The formulas analogous to Eqs. (|5.ip and (|5.4p 
for the projected distributions, denoted 'hQ(x*,t*) and 
fiMix* ,t*), are 



V. SPATIAL DISTRIBUTION FUNCTIONS 

Yet another view of the normal and anomalous dif- 
fusive behaviors is obtained by looking at the spatial 
distribution functions themselves, that is, by comput- 
ing nG{r*,t*) and nj\/(r*,t*). The Fourier transforms in 
Eqs. (|3.9p and p.lOp can be inverted analytically, yielding 

nGir*.n = Sir*) " ^ ^ rdr^t*), (5.1) 

where 



2tt iu 



and 



(5.2) 



k{u) ^ \ 2 l-iPaiu) +u/A , ReK>0. (5.3) 



nG{x*,t*) = cbcit*) 5{x*) + - TG{x\t*), (5.6) 



and 



(5.7) 



where the functions Fq and Fm are the same as those 
given in Eqs. (|5.2p and l|5.5p with r* replaced by x* . 

The 5 functions in ^(^(r*,^*) and nG{x* ^t*) are the 
results of our having neglected the intra-cage motion dis- 
cussed earlier. In this approximation, with spatial res- 
olution only of the order of £, molecules in frozen do- 
mains remain exactly at their initial positions with de- 
caying probability (j)G(t*). The distributions nMir* ,t*) 
and nM{x*,t*) have no (5-function contributions because 
diffusion starts immediately in mobile regions. In graphs 
of experimental or computational data, these S functions 
are visible as narrow Gaussian displacement distributions 
centered at a;* = 0, whose root-mean-square widths, say 
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l'^ = ic/R*, are the average intra-cage displacements 
that we are neglecting here. 

As a first look at the spatial distribution functions, 
in FiglU we show log^Q[Fs{x* ,t*)] as a function of the 
displacement x* for a sequence of times t* . In analogy to 
Eq. ()4.1|) . we write 

Fsix*,t*) ^VGnGix*,t*) + il-VG)nMix*,t*), (5.8) 

For plotting these graphs, and for this purpose only, we 
replace the d function in Eq. (j5.6p by a normalized, Gaus- 
sian, intra-cage distribution 



S{x*)^pcix*) 



1 



(27r^^2)i/2 



exp 



2£% 



(5.9) 



We make this replacement primarily because the result- 
ing graphs look - and indeed are - more realistic this 
way; but we emphasize that this is not a systematic cor- 
rection of the small-time behavior because the intra-cage 
fluctuations are not otherwise included in nQ{x*,t*) or 
nM{x*,t*). Again, we choose A = 100 and Vg = 0.5; 
and, so that the small- a;*, early-time behavior be visible 
in the figure, a relatively large but possibly realistic value 
of = 0.1. Because of the factor (pG{t*) in Eq.(|5.6p. this 
central peak in hG{x*,t*) disappears after times of the 
order of r*. 

At small times t* , and for sufficiently large displace- 
ments cc*, the graphs in Figd] exhibit the broad exponen- 
tial tails reported elsewhere in the literature (see (3l.[l3|l. 
The exponential tail is a robust mathematical feature of 
this class of models, closely associated with the decou- 
pling of the mobile and glassy behaviors, and not depen- 
dent on any details of the glassy waiting-time distribution 
^pG{t*)- The degree of decoupling is reflected here by the 
magnitude of the parameter A, which must be large for 
strong decoupling. The only requirement on tpGit*) is 
that it decays rapidly enough that ^g{u) becomes van- 
ishingly small at large u. 

To see how the exponential function emerges, evaluate 
the M- integrations in Eqs. (|5.2p and (|5.5p by first integrat- 
ing around a circle of radius uq centered at the origin, 
and then closing the contour around the cut on the neg- 
ative M-axis from u = —uq to u —oo. (The curves 
in Fig|4] were computed numerically with uq = 2.) If 
A ^ 1, we can choose 1 ^ uq <C A, so that i^G{u) ~ 
and k{u) « everywhere around the circle. Then, 
for small i*, the integration around the circle yields 
TG{x*,t*) w TM{x*,t*) w exp(-V2a;*). This limiting 
behavior is most apparent in Figd] for the two earliest 
times, t* = 0.03 and 0.1, where the slope has the pre- 
dicted value of — \/2. 

At larger times, and in the limit of indefinitely large A, 
the preceding argument implies that all of these curves 
approach the same slope at large x* . On the other hand, 
at very early times t* < A~^, the integration is domi- 
nated by the discontinuity across the cut near u = —A, 
and the displacement distribution is dominated by the 
early diffusive motion of the initially mobile molecules. 



Another interesting case is the limit of large t* and 
X* <^ t*. The integrands in Eqs. (|5.2p and (|5.5p vanish 
like exp (— l/w) in the limit —u = w — > 0+. Therefore, we 
can use the long-time, small-w approximation, A{w) « 
1 + 3w/2, to estimate 



Tg{x* ,t*) K Tm{x* ,t*) (X e^p 



3x* 



At* 



(5.10) 



thus the spatial distribution reverts to a Gaussian. 
Ea. (|5.10p fits the curve in Fig|4] for t* = 10 reasonably 
well for x* < 5. 



VI. MOMENTS OF THE SPATIAL 
DISTRIBUTIONS 

Another way of looking at these spatial distributions 
is to compute their time-dependent moments (r*'"(i*)). 
We use the identities 



and 



{r*\t*)) = 15 



dk'^ 



k=0 



k=Q 



(6.1) 



(6.2) 



When inverting the Laplace transforms as in Eq. 
we again close the contour of integration in the negative 
u plane, but in this case it is mathematically essential to 
include the circle around the origin in order to account 
for the singularities that occur there. The mathematical 
situation becomes clear by inspection of the formula for 



{r*\t*))M - 6 



du 



2m k'^{u) 
du 



2m 1 _ iI}g{u) + w/A 



(6.3) 



Because ipGiu) ~ 1 — 3u/2 for small u, this integrand has 
a first-order pole at the origin. Higher-order moments 
have higher-order poles. 

We can use the analogs of Eq. (|6.3p for higher moments 
and for both the ug and hm distributions to obtain ex- 
act asymptotic results in the limits of very small and 
very large t*. For vanishingly small values of t* , close 
the contour around a circle at large values of u where 
i/jGiu) ~ 2/u. In this limit, the integrand behaves like 
exp(Mt*)/M^, and 



(r* 



3 At* 



Similar small-t* calculations yield: 



(r**)M« 15 A^r^ 



(6.4) 



(6.5) 
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FIG. 5: (Color online) logiQ[(r* {t*))F] as a function of 
logio(t*) for Vg = 0, 0.5, 0.75, 0.9, 0.97, 0.999, and 1.0 (red 
curve), reading from top to bottom. 



and 



3 At* 



10 t*^. 



(6.6) 



(6.7) 



For very large t* , close the contour in Ea. (j6.3p on 
a vanishingly small circle around the origin and use 
tpciu) ~ 1 — 3 u/2. The result is that, after times so long 
that the tagged molecule no longer remembers whether it 



started in a mobile or a glassy region, both 



/M 



and 



(r* )c converge to the same slowly diffusing Gaussian 
distribution for which (r*^) « 2t* (with corrections of 
the order of 1/A). 

Note several features of these results. As expected, the 
initially mobile molecules exhibit rapid (large A) Gaus- 
sian diffusion at small times. On the other hand, the 
displacements of the initially frozen molecules are non- 
Gaussian and "pseudo-ballistic" with (r*^)^ ~ At*'^. 
This behavior has nothing to do with early-stage bal- 
listic motion of molecules within their cages but, rather, 
is an intrinsic feature of intermediate-stage, anomalous 
diffusion in this model. 

We turn finally to the full time dependence of the mo- 
ments. The results for the weighted average, 

{r*\t*))F = Vg {r*\t*))G + {l-VG) {r*\t*))M. (6.8) 

are shown in Figl5]for a series of different values of Vg^ 
as functions of log]^Q(i*). At very short times, these mo- 
ments rise as functions of t* according to the estimates 
in Eas. (|6.4p - (|6.7p . The pseudo-ballistic behavior is ap- 
parent only for Vg ~ 1 because, for smaller Vg, these 
weighted moments are dominated by the displacements 
of the initially mobile molecules. After the initial rise, at 
roughly t* « A~^ = 10~^, the moments cross over from 



FIG. 6: Non-Gaussian parameter 02 as a function of logig(t* 
for Pg = 0. 



fast (3 relaxation to slow a decay as the tagged molecules 
are repeatedly trapped and then escape from glassy do- 
mains. Finally, in Fig|Sl we show the non-Gaussian pa- 
rameter 



a2{t*) 



M 



'•2(t*))2 



1 



(6.9) 



M 



only for Vq = because, in this approximation where 
we have neglected intra-cage vibrational motions, only 
the mobile molecules exhibit Gaussian displacement dis- 
tributions at early times. Here we see explicitly that the 
non-Gaussian behavior occurs during the crossover from 
/3 to a relaxation. Both of these last two figures are ar- 
tificial in the sense that we are independently varying 
the parameter Vg, which in fact ought to be determined 
uniquely by the temperature along with the parameters 
R* and A. 



VII. CONCLUDING REMARKS 

The model proposed here contains a relatively simple 
mechanism for producing stretched exponential decay of 
molecular correlations, as observed via the self interme- 
diate scattering function. That mechanism emerges di- 
rectly from the spatial heterogeneity of a glass-forming 
liquid. The model illustrates how anomalous diffusion, 
as exemplified by a broad exponential (non-Gaussian) 
tail of the molecular displacement distribution, is related 
to heterogeneity and ~ indirectly - to the accompanying 
stretched exponential behavior. Yet another feature of 
the model is the way in which it illustrates how the tran- 
sition from relatively fast (3 relaxation to slow a decay 
correlates with the onset of anomalous diffusion. Our 
analysis neglects the very short time, strongly localized, 
intra-cage fluctuations. The transition between the latter 
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motions and cage-breaking events seems to us to be out- Acknowledgments 
side the range of validity of the continuous-time random- 
walk theory used here. 

There are, of course, many missing ingredients, almost 
all of which pertain to the temperature dependence of 

the various parameters introduced here. One of us (JSL) This research was supported by U.S. Department of 

plans to address those issues in a subsequent report. Energy Grant No. DE-FG03-99ER45762. 
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